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Abstract 

We apply QCD Finite Energy Sum Rules to the scalar-isoscalar current to deter- 
mine the lightest uu + dd meson in this channel. We use 'pinch-weights' to improve 
the reliability of the QCD predictions and reduce the sensitivity to the cut-off sq. 
A decaying exponential is included in the weight function to allow us to focus on 
the contribution from low mass states to the phenomenological integral. On the 
theoretical side we include OPE contributions up to dimension six and a contribu- 
tion due to instantons taken from the Instanton Liquid Model. Phenomenologically, 
we incorporate experimental data by using a coupling scheme for the scalar current 
which links the vacuum polarisation to the vrvr scattering amplitude via the scalar 
form factor. We find that the sum rules are well saturated for certain instanton 
parameters. We conclude that the /o(400 — 1200) definitely contains a large uu + dd 
component, whereas the /o(980) most likely does not. We are able to estimate the 
average light quark mass and find mq{l GeV^) = 5.2 it 0.6 MeV. 
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1 Introduction 



Currently, experiment suggests that there are more scalar mesons than can be accommo- 
dated in a single quark-model nonet. Many QCD motivated models predict the existence 
of non-gg mesons such as glueballs hybrids 0, qqqq [|| and KK states. It is 
tempting to say that these 'spare' experimental states are the unconventional mesons al- 
lowed by QCD but, with scalar glueball and multiquark states expected to have masses 
comparable to the conventional qq states, which are the extra states? In the work that 
follows we hope to shed some light on this question by attempting to determine which of 
the 5 scalar-isoscalar mesons currently listed in the Particle Data Group tables is the 
lightest uu + dd meson in this channel. To address this question we will use the QCD 
sum rule technique. The answer has importance for chiral symmetry breaking in QCD, 
see for instance |^. 

QCD sum rules are integral expressions that relate the hadronic and partonic regimes, 
i.e. the low energy world of resonances with the high energy world of (tractable) QCD. 
Since their conception over twenty years ago |p, QCD sum rules have become an estab- 
lished technique both for calculating hadronic properties in channels, where the QCD 
expressions are under control, and, conversely, estimating QCD parameters (such as the 
masses of the quarks) in channels where there is good experimental information. 

Using QCD sum rules in the scalar channel is not a new idea. They were first applied 
to the scalar mesons in |^ where Laplace sum rules were used (for a recent review of 
Laplace Sum Rules see llO|)- Here the phenomenological side was represented using 
the, now standard, resonance + continuum ansatz, with resonances being represented as 
(5-functions and the continuum being calculated entirely within perturbative QCD. As 
only the Operator Product Expansion (OPE) terms, discussed in Sect. |, were taken 
into account on the theory side, exact mass degeneracy between the hghtest isoscalar 
and isovector was found, with rrifg = rriag = 1.00 ± 0.03 GeV. The ss state was predicted 
to have a mass of around L35 GeV. These findings were supported by Bramon and 
Narison [|n|], who used QCD sum rules to calculate the couplings of the /o(980) and 
ao(980) to two photons. Again modelling the resonances with 5-functions the authors 
concluded that a qq interpretation of these two mesons could not be ruled out with the 
data then available. 

The effects of going beyond the OPE were studied in [jT^, where instantons (see 
Sect. ^) were included on the theoretical side, but not in the QCD continuum on the 
phenomenological side. Once again the resonances were treated as 5-functions and only 
the isospin-1 channel was considered. It was concluded that the lightest qq state in this 
channel had a mass < 1 GeV. 
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For the light scalar mesons the zero-width approximation is not a good one and the first 
attempt to go beyond it in a sum rule investigation was by Elias et al |TB|. Here resonances 
were represented by Breit-Wigner formulae, which, whilst better than a 5-function, still do 
not adequately describe the complex structure in the scalar sector. A further abstraction 
was introduced by replacing these Breit-Wigner shapes with a Riemann sum of rectangular 
pulses. Laplace sum rules were used, but in an integral form which requires the calculation 
of perturbative expressions at rather low energy scales. Again instantons were included, 
but not in the continuum. In the isovector channel, the ao(980) was found to decouple 
from the sum rule and the mass of the lightest quarkonium isovector was consistent 
with the ao(1450). In the isoscalar channel the conclusions were not so clear cut, it was 
predicted that the lightest quarkonium here should have a width less than half of its mass, 
and although the /o(980) could not be ruled out as the main contributor a lighter state 
was preferred. Similar conclusions were reached in This study used both Laplace 
and Finite Energy sum rules, again requiring perturbative expressions to be evaluated 
at low energies. The authors noted, that for a consistent treatment of the Borel sum 
rules, the instanton effects should also be included in the QCD continuum contribution. 
The sum rules were dominated by an isoscalar with mass around 1 GeV and an isovector 
with mass around 1.5 GeV, thus suggesting the /o(980) and ao(1450) as the lightest qq 
states in their respective channels. However, in the same work, a comparison to a more 
realistic resonance shape predicted resonance parameters for the lightest quarkonium of 
m ~ 860 MeV and F ~ 340 MeV. These parameters are not consistent with the /o(980), 
but could describe a Breit-Wigner fit to the /o(400 — 1200). 



In [0, Maltman used pinched weight finite energy sum rules to calculate the decay 
constants of the ao(980) and ao(1450), which were then compared with the decay constant 
of the, presumably qq, ii'g(1430). Instantons were included and all QCD integrals were 
carried out in the complex plane, which improves the convergence of perturbative expres- 
sions. The isovector decay constants were found to be of comparable size, suggesting a 
similar structure for both. The author concluded that this favoured a Unitarised Quark 
Model [0 scenario, where two physical hadrons can arise from one 'seed' state. 

In the work that follows we will apply QCD Finite Energy Sum Rules introduced in 
Sect. 2 to the scalar-isoscalar channel. On the theoretical side we will include in Sects. 3 
and 4 instanton effects as well as the OPE terms up to dimension six and use the contour 
improvement prescription to increase the convergence of perturbative expressions. On the 
the phenomeno logical side in Sect. 5 we will incorporate experimental data directly by 
relating the hadronic vacuum polarisation to the tttt scattering amplitude via the pion's 
scalar form-factor. The sum rules are evaluated in Sect. 6 and we present our conclusions 
in Sect. 7. 
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2 The Sum Rules 



The connection between quarks and hadrons is through the vacuum polarisation, n(s). 
This is defined by the correlation function of two currents, i.e. 

nw^./d'..-(0|T{.W.(0)}|0). (1, 

where s = = —Q^ and the currents are chosen so as to select the desired channel. In this 
work we are considering the scalar-isoscalar mesons and so we choose the Renormalisation 
Group invariant current 

J{x) = — ^ {u{x)u{x) + d{x)d{x)^ = -—^n(x)n{x) , (2) 
2 v2 

where rrig = |(m„ + m^) and we use the symbol n to denote an effective light-quark. 
The second equality follows because we ignore the two-gluon intermediate state, whose 
contribution is 0{rn?^) and so can be safely neglected. 




s 



Figure 1: The complex s-plane showing the so called 'Pacman' contour which is 
used to derive the FESRs. 



n(s) is analytic everywhere in the complex s— plane except along the time-like axis, so 
for any weight function, w{s), which is real for all real s and analytic within the 'Pacman' 
contour shown in Fig. |I| we can write 

- f w{s) Imn(s) ds = — ^— l w{s) n(s) ds . (3) 

On the left-hand side of Eq. we have an integral in the physical region, and here we 
choose to describe the integrand using hadronic physics (either with experimental data or 



4 



a phenomenological parameterisation). On the right-hand side we have an integral in the 
non-physical region (complex s), and here we use field theory [i.e. QCD) to calculate the 
integrand. The assumption that Eq. (|^) holds for a range of values of sq is a statement 
of global quark-hadron duality, i. e. the integral of an hadronic expression is equal to the 
integral of an expression written in terms of partonic degrees of freedom. 

Early sum rule work used positive integer powers of s as weight functions. In this work 
we wish to investigate the lightest meson in the I = J = channel, and so we choose a 
decaying exponential, w{s) = exp(— s/M^). On the phenomenological side, this weight 
function will suppress the region of the integral s > M^, and so by varying we can 
control which regions of the integrand contribute most significantly to the sum rule. If 
we also include integer powers of s in the weight function then we can build up a family 
of sum rules. 

Almost by definition, at energies where resonances dominate the spectral density, the 
implications of QCD are not straightforwardly calculable. If QCD were tractable then 
the technology of QCD sum rules would be redundant. This means that, for moderate 
values of sq, our straightforward QCD description of n(s) must fail, at the very least, 
in the region of the positive real axis (see Fig. ||). However, it is believed [17| that it is 
only in the region of the positive, real axis that the straightforward QCD description is 
particularly unreliable and that in the rest of the s-plane calculable QCD provides a good 
description of the correlator. Consequently we introduce a zero into the weight function at 
the point where the circle joins the real axis, i.e. Sq. Thus our sum rules have the weight 
function w{s) = {1 — s/sq) e~'^/*^^ and we denote them by the symbol T^. Elsewhere in 
the literature they are often referred to as 'pinched weight' sum rules. Then the left hand 
side of Eq. (3) becomes 

Tfc(so, M^) = ^j\s^ (^1 - exp (-^) Imn(s) ds . (4) 

Introducing this zero into the weight function has a number of advantages. Firstly, on the 
theoretical side we avoid the problem, just described, of the failure of QCD on the part 
of the circle near the positive real axis. Secondly, on the phenomenological side, we have 
introduced a second factor with a tendency to suppress the higher energy portion of the 
integral and consequently reduce the dependence on the exact value of Sq [T^, [T^ . 

The parameter sq is not entirely free: it must be large enough that the QCD expressions 
we are able to calculate are expected to be a good approximation to the full correlator, yet 
So should be small enough that we have experimental information. The art of sum-ruling 
is to find (sensible) values of sq for which [Tfc(so, M'^)]had ~ [^^(-50, ^'^)]qcd ^^^^ ^ wide 
range of the unphysical parameter M^. When this occurs we say that the sum rules are 
saturated. 
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Although the current Eq. is RG- invariant, the correlator it gives rise to is not. As 
well as making the calculation of the correlator within QCD more difficult, this would also 
lead to an unwanted scale dependence in our final results. The second derivative of this 
correlator, Il"{s), is RG-invariant. As n"(s) is also analytic within the contour shown in 
Fig. |1], we could choose to write the sum rules completely in terms of this quantity, but 
Imn"(s) is not directly related to any physical quantity and so we would lose our link 
to experiment. Instead we use partial integration to re-express the theoretical side of the 
sum rules in terms of n"(s). To do this we write the integral on the theoretical side as 

ii = ^v"{s)U{s) ds . (5) 

We now integrate this expression by parts twice to give, 

ii = v'{s)U{s) - v{s)n\s) + j v{s)Tl"{s)As , (6) 

and arrange the constants of integration which appear in v{s) and v'{s) so that these 
functions disappear at s = sq. Then our sum rule is 

i w{s) [Imn(.)],^, d. = £_^i;(.) [n"(.)]Q^^ d. . (7) 



3 The Operator Product Expansion 

In order to calculate the vacuum polarisation within QCD we make use of the Operator 
Product Expansion (OPE) ||20[. This replaces the matrix element of the product of 
currents in Eq. (|1|) by a linear combination of local operators ordered in terms of increasing 
dimensionality, 

Tl{s=-Q') = Y,Cn{Q'){On) . (8) 

n 

The first term in this expansion is the purely perturbative contribution, with higher 
dimension terms giving the contribution due to various condensates. We will consider 
terms up to dimension six, higher order terms are expected to be negligible. Expressions 
are available in the literature for the dimension zero term to four-loop order, the dimension 
four term to two loops and the dimension six term at leading order. We quote these 
expressions with refernces to the original works for the full information. All expressions 
are given in the MS scheme. 
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The dimension zero expression can be found in [p^ and is given by 
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and ({n) is the Riemann Zeta function. As n"(s) is renormalisation group invariant, we 
are free to make any convenient choice of the renormalisation scale. We use the 'contour 
improvement' prescription, whereby the logs in Eq. (|^) are 're-summed' by making the 
choice = Q"^. This improves the convergence of the perturbative series by removing the 
dependence on large powers of possibly large logs which would appear in the (unknown) 
higher order terms. The contour improved perturbative contribution to our correlator is 
then, 

'5071 
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144 



fc(3, 



9 6 4 ^ 



(10) 



Naively, the only gauge-invariant operator of dimension two that can be constructed 
within QCD is the so-called mass insertion term, m^. For the light quarks that we are 
considering this term can safely be ignored. Some recent works have suggested that 
there may be mechanisms outside of the normal OPE, which could lead to dimension two 
contributions to correlation functions. These mechanisms include renormalons and 
a tachyonic {i.e. imaginary) gluon mass p3[. It has been argued that, if they do exist. 



the effects of dimension two operators would not be significant in most sum rule analyses 
and that even in calculations where they would be expected to be most noticeable, these 
effects are consistent with zero [24]. Hence, in this work, we choose to ignore them. 



The contour-improved contribution at dimension four can be inferred from [^, ^ ^ 
and is given by 

22a,(Q2 



2Q6 



6 {nigqq) 
4a,(Q 



1 + 



+ 1 



9 



3 

287r2 



9 



1 + 



121a,(Q2 



(11) 



7 



where 



= 1 + ^ +4a,{m,ss) l + ^ l + ^ , (12) 



24 / 47r2 V 3 
and 

The leading order dimension six contribution can be taken from and we follow this 
reference in taking this term at the fixed scale of 1 GeV. The dimension six contribution 
is then 

„ , 3m2(l GeV^) / 



+ TT^a 



Q 

2 
3 



(14) 



where 7^ are the Dirac matrices, a^^ = |{7^, Ju} and A" are the usual Gell-Mann matrices 
normalised such that Tr[A'*A^] = 25"*. Currently there is no reliable method to calculate 
the vacuum expectation values appearing in Eq. ([T^) , so we follow standard practice 
and relate these dimension six condensates to the (RG variant) dimension three light 
quark condensate (qq). For the mixed condensate, {gqa^uG'^'^"'X"'q) , this is done via the 
parameterisation 

(gqa^M^^'X^q) = rnUqq) , (15) 

and for the four-quark condensates, {qa^i,X"'qqcr^'^ X'^q) and {q-yfj_X'^qqY'X"'q), via the vac- 
uum saturation hypothesis [||. To take into account the possible violation of the vacuum 
saturation hypothesis we include a multiplicative factor 14s- Thus our final expression for 
the dimension six contribution is 

nW^)^g(»SK5,>-^^^(5,>^) . (16) 

In order to evaluate all these expressions we need to calculate the running coupling and 
quark masses at complex values. This is done by solving the four-loop Renormalisation 
Group Equations that define the QCD (3 and 7 functions (our conventions for these 
functions are given in Appendix |^). 

For the coupling, we choose to solve the implicit equation 

da' , f 
^7ZK - log ( : 

as(m2) 



f3{a') \m: 



log ( ) = , (17) 
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iteratively using Newton's method. As input, we take the experimental value as{iTL^) = 
0.334 ± 0.022, as measured by the ALEPH collaboration For reference, had we 

used the standard expansion for as, then at four-loops and with three active flavours, 
this central input value would correspond to Aqcd ~ 365 MeV. For higher values of 
as(m^), or correspondingly Aqcd, the radiative corrections to Eq. ( ]T0| ) become more 
important and start to swamp the higher dimension terms, i. e it becomes more important 
to include higher loop corrections to the perturbative contribution than to include the 
condensates. It has been argued that the sum rule methodology could become invalid 
for Aqcd ^ 330 MeV. This is especially problematic for attempts to determine condensate 
values via the sum rule technique. As this is not our aim and considering the success of 
previous sum rule studies with similar and higher values of a;s(m^) (and Aqcd), we feel 
justified in applying sum rules to this problem. 

The momentum dependence of the quark mass is given by 



mq{jj?) = mg{l GeV^) exp 



a^lGeV^) 



7(aO 



da' 



4 Instantons 



The OPEs for the scalar-isoscalar and scalar-isovector channels are identical. This was 
originally interpreted as an explanation for the almost exact mass degeneracy of the 
/o(980) and the ao(980) However, we know that the OPE is not enough to describe 
the QCD vacuum completely. Other effects are important, particularly in the scalar and 
pseudoscalar channels |3^ . 

The complex nature of the QCD vacuum means that it is not currently possible to solve 
the equations of motion of the full theory directly. We must make use of approximations, 
which in this work we choose to be the Instanton Liquid Model Here the 

QCD vacuum is modelled as a four- dimensional 'liquid' of instantons, which is assumed 
to be dilute enough to make the idea of individual instantons meaningful. The liquid is 
characterised by the average instanton size pc and four- volume density ric (or alternatively 
the average instanton separation, Vc = Uc ^^^). The Instanton Liquid Model is then thought 
to be a good approximation, if the ratio pc/rc is small. Within this model the instanton 
contribution to the scalar-isoscalar light-quark correlator was found in [^, |3^ to be 



27r2 



i^i(PcV^ 



(19) 



where Ki{x) is a modified Bessel function of the second kind (or MacDonald function). 
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Taking two derivatives we obtain 



nr(.) 



3pg ml 
47r2 



(20) 



As stated earlier n"(s) is an RG-invariant quantity. The only quantity in Eq. (pO|) that 
could have a renormalisation scale dependence is the quark mass. This renormalisation 
scale dependence cannot be cancelled out anywhere else, so we conclude that the quark 
mass appearing in Eq. ( ]20|) must be the quark mass at some fixed scale, v"^. A priori we 
do not know what this scale is. It is reasonable to expect that it is related to the instanton 
scale given by 1/pc oi the hadronic scale of ~ 1 GeV, but we have no physical reason for 
choosing 1/pc rather than e.g. Ijpc- Thus we take this scale as an input parameter in 
the calculation. If this scale is found to be considerably different from the hadronic and 
instanton scales, then it would call into doubt our treatment of the instanton contribution. 



5 The Phenomenological Side 
5.1 The Coupling Schemes 

At low energies the only hadronic process possible in the scalar-isoscalar channel is tttt 
scattering. As pions are made of non-strange quarks, it is reasonable to expect that, 
in the elastic region, the spectral density for the correlator is related to the absorptive 
part of the vrvr tttt scattering amplitude, T{s). This link can be made via the scalar 
form- factor of the pion, d{s), which is defined by 

d{s) = (0|m,gg|7r7r) . (21) 

Inserting a complete set of hadronic states into the correlator Eq. (|l|) and keeping only 
the lowest state, i.e. two pions, gives the spectral density in the elastic region, 

Imn(s) 3/3(s), ,,2 



TT 327r2 



\d{s)\' , (22) 



where P{s) = yl — Aml/s. 

Watson's Final State Interaction Theorem |^ implies that the form factor and the 



scattering amplitude have the same phase. Thus, in the region where elastic unitarity 
holds, we can write 

T{s) = a{s)d{s) , (23) 
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where the couphng function, a(s), is real for real s > 4m^. Substituting Eq. ( p3D 
into Eq. ( p2|) and making use of the elastic unitarity relationship ImT = /3 |Tp, we obtain 

Imn(s) 3 ImT 

It is well known that vrvr scattering contains an Adler zero , and that this zero does 
not appear in the scalar pion form factor. We now make the simplifying assumption that 
this is the only s-dependence appearing in a{s), i.e. 

a{s) = ao{s - sa) (25) 

where sa is the position of the Adler zero and ao is a constant. Substituting Eq. (^) into 
Eq. (p^) gives our first ansatz 

Imn(g) ImT(g) 

= / 7 , (26) 

TT (s - Sa)^ 

where / is an unknown constant of proportionality. 

Although the Air threshold is at 560 MeV, it is well known that multi-pion channels 
are not significant below around 1400 MeV, when quasi two-body channels, e.g. aa, pp, 



become important ||3^, |3^. The first important inelastic channel is vrvr KK. Even after 
this has opened the vrvr final state will continue to pick out the non-strange contribution 
to the sum over states. Thus we might expect Eq. (|26| ) to give a reasonable approximation 
to the true spectral density even above 1 GeV. We call this Coupling Scheme I. 

The multi-pion final states mentioned above will increase the nn spectral density. 



causing Eq. (^) to be an underestimate at higher energies. Phenomenologically, we 
might expect the true spectral density to be enhanced by the ratio of the total to elastic 
vrvr cross-sections, thus giving 

which we call Coupling Scheme II. In the elastic region, the two coupling schemes are, 
of course, identical. Above the KK threshold, whereas Coupling Scheme I ignores all 
inelastic channels and so will underestimate the spectral density. Coupling Scheme II 
takes into account all inelastic channels, including those with hidden strangeness, and so 
may be an overestimate. 

As input we take the parameterisation of the / = J = vrvr scattering data carried out 



by Bugg et al ||39|. However, when continued down to threshold this fails to reproduce 



the scattering length predicted by Chiral Perturbation Theory (xPT) and so in the region 
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Figure 2: A sketch showing the absorptive part of the / = J = tttt scattering 
amphtude as calculated from the data set described in the text (line). Superim- 
posed are experimental points from [^] (diamonds), (boxes), (triangles) 
and E5| (circles). 



below 450 MeV we use the Roy equation extrapolation of the classic Ochs- Wagner phase- 
shifts [^, as carried out by Pennington and used in the discussion of 77 tttt ||41|| . 

Fig. ^ shows a sketch of the absorptive part of the isospin zero, S-wave amplitude 
for TTTT scattering as calculated from our data set. From this we can clearly see that the 
/o(400 - 1200) and the /o(1370) appear as peaks, while the /o(980) and the /o(1500) 
show up as sharp dips. Thus these latter two states, and in particular the /o(980), will 
largely decouple from the spectral density and if saturation of the sum rules is achieved 
with either of our two Coupling Schemes, it will not be due in any great part to the 
/o(980) or /o(1500). By varying both and Sq it may be possible to determine which 
of the /o(400 — 1200) and the /o(1370) plays the dominant role in saturation, but for the 
evaluation of the Wilson coefficients in the OPE to be trustworthy, we would expect sq 
to be too high to exclude the /o(1370). Nevertheless, as the decaying exponential in the 
weight functions of our sum rules suppresses the integrand for s > M^, it is still possible 
to determine which resonance is most important for saturation. 
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5.2 Normalisation 



A value for the normalisation factor /, can be determined by making a comparison with 
the spectral density obtained using an Omnes representation |^ of the form factor. This 
representation is given by 

(h(s') ds' 



d{s) = d{0) exp 



vr J4^2 s'{s' 



i2i 



where is the phase of / = J = tttt scattering. The form factor at zero momentum 
transfer, which normalises the Omnes function, can be determined from the Feynman- 



Hellman Theorem and is given by 

rf(0) 



dnin 



(29) 



To leading order in xPT, the physical pion mass is given by the Gell-Mann-Oakes-Renner 
relation, A{mqqq) = —m^fl [|8|, and thus we see that 



d{Q) 



(30) 
01) of the 



Going to next-to-leading order in xPT would introduce a correction to Eq 
order of 1% (see for example ^^), which is safe to neglect. 

To evaluate the integral in Eq. (|28| ) exactly requires knowledge of the scattering phase 
in the high energy region where it is unknown. However, the integrand is dominated by 
the region around s. So if we restrict our comparison to the low energy region then the 
behaviour of the phase above some point Si, where si ^ s, will not greatly affect the 
value of the integral. We take Si to correspond to the last data point and approximate 
the phase above this point as a constant, (pi. Thus, the spectral density, as calculated 
from the Omnes representation, is written 



Imnfs) 



TT 



3m: 
32^ 



Si — s 



20i/7r 



exp 



2s 

TT 



0(s')ds' 

4m2 S'{s' — S)_ 



(31) 



where the bar on the integral sign denotes the Cauchy Principal Value. 

The value of the constant, /, can then be determined by fitting Eq. (|26D to Eq. (pT|). 
The position of the Adler zero for isospin-0, S-wave vrvr scattering is fixed to be close to 
which is its position at lowest order in xPT and so we essentially only have one 
parameter, /, to determine. The fit, shown in Fig. |^, was carried out up to a maximum 
energy of 400 MeV. Over this energy range it does not matter which Coupling Scheme is 
used. Each gives the value / = 6.9 x 10^^ GeV^. However, this value is quite sensitive 
to the maximum energy of the fit and should not be trusted to an accuracy better than 
20%. 
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Figure 3: The spectral density below 400 MeV as calculated from our Coupling 
Scheme (points) and the single channel Omnes representation (solid line). 



6 Results 

In any sum rule calculation, we hope to find regions of the parameter space where the 
theoretical and phenomenological sides are equal. We measure this saturation of the sum 
rules by introducing the double ratios 

D.iso, M^) = ^^±^1^^ . (32) 

This ratio is independent of the overall normalising constants on either side, i.e f and 
mg(l GeV^), and obviously is equal to one when the sum rules are saturated. As sum rule 
analyses are generally accurate to around 10-20% we set an agreement region of (1 ±0.1). 
In practice we shall use only the lowest such double ratio for each type of sum rule, i.e Dq, 
as the integer powers of s in the higher sum rules increase the importance of the upper 
end of the integrand on the phenomenological side. 

Our sum rules contain a number of input parameters which we list in Table along 



with the values used. Any deviations from these values will be discussed in Sects. 5.1 



^ and |6.3|. Where appropriate we also include a brief description of where these values 



come from. 
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Parameter 


Value 






0.334 


From Experiment |^ 


{mg-qq) 


-(95MeV)4 
(381 MeV)^ 
-(195 MeV)^ 


From GMOR relation [|§ 






{ttIsSs) 




m,(l GeV^) 


159 MeV 






-(225 MeV)3 
0.8 GeV^ 






Sum rule Estimate ^0 




1 


Pc 


(600 MeV)-i 




V 


~ IGeV 


To be determined 



Table 1: Standard values of all parameters used to evaluate the sum rules. 
6.1 The Instanton Parameters 

The first task is to determine ^ the scale of the quark mass in the instanton contribution. 
To evaluate Eq. (^), we introduce a new parameter, A, which is equal to the ratio of 
m2(i/2) to m\{\ GeV^), i.e., 



A 



mg(z/2) 



mg(l GeV 



exp 



as(lGeV2) 



7(aO 
/3(a0 



da' 



(33) 



where the second equality follows from Eq. ([T8|) . In Fig. ^ we show the stability curves 
obtained for various values of this parameter using Coupling Scheme II and sq = 3.7 GeV^. 
We have found that the saturation of the sum rules is practically independent of Sq in the 
region 3 GeV^ < So < 4 GeV^ and so no significance should be attached to this choice of 
Sq. However, it was found that saturation could not be achieved with the standard value 
of Pc = (600 MeV)~^, and in Fig. | we pre-empt the results of Figs. ^ and |^ by using the 
non-standard value pc = (480 MeV)^^. 

It is immediately clear that the sum rules are far from satisfied when A = 0, i.e. when 
the instanton contribution is removed. However, the sum rules are well saturated for 
A ~ 1.2 and we have found that the best value is A = 1.25, which we will use from now 
on. This value of A corresponds to a scale u = 0.96 GeV. 

As stated earlier saturation can only be achieved by use of a non standard value of 
Pc- In Figs. ^ and ^ we show the effect that varying this parameter has on the saturation 
curves. We see that good saturation is achieved for 1/pc = 470 — 490 MeV, with the 'best' 
value being pc = (480 MeV)^^ ~ 0.42 fm. This value is slightly larger (25%) than the 
estimate |^ based on the size of the gluon condensate which has become the standard 
value employed in the Instanton Liquid Model. However, it is entirely consistent with 
current lattice determinations which tend to find sizes in the range 0.32 - 0.43 fm and, 
unless stated otherwise, we will use this 'best' value of pc from now on. 
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0.2 - 

I I I I I I I I 

0.5 1 1.5 2 2.5 3 3.5 4 

M2 (GeV2) 

Figure 4: The saturation curves obtained for various values of A, with = 
(480 MeV)~^ and all other parameters as listed in Table 0. 

6.2 Saturation 

As shown in Fig. |], once the instanton parameters have been determined the quality of 
the saturation of our sum rules is very good. This means that the /o(980), being largely 
decoupled from our Coupling Schemes, does not play an important role in this saturation. 
On the other hand, the saturation curves are seen to be flat from quite low values of M^. 
With our weight function suppressing the portion of the integrand s > relative to the 
low energy region, the shape of the curves implies that saturation is achieved long before 
the /o(1370) contributes significantly to the phenomenological integral. This suggests 
that it is the /o(400 — 1200) that plays the most important role in saturating our sum 
rules. 

As the values of the condensates have not been fixed definitively, we have performed a 
number of tests to determine how sensitive the saturation of our sum rules is to variations 
of these input values. To test the dependence on the size of the gluon condensate, we 
write 

(asG^) = c{asG^)o (34) 

where (asG'^)o is the standard value listed in Table |l|. In Fig. |^ we show the saturation 
curves for various values of this factor, c. We see that the saturation is not strongly 
dependent on the size of the gluon condensate, but the best results are obtained in the 
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1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 




i' \ ' 


1 1 1 1 1 

= (400 MeV)-^ 




= (450 MeV)-^ 




p^ = (500 MeV)-^ 




p^ = (550 MeV)-^ 




p^ = (600 MeV)-^ 





0.5 



1.5 2 2.5 
M2 (GeV2) 



3.5 



Figure 5: The saturation curves obtained for various values of pc, with A = 1.25 
and all other parameters as listed in Table |I|. 



p. = (470 MeV)- 
p^ = (480 MeV)-i 
pZ = (490 MeV)-^ 




Figure 6: As Fig. ^ but for finer steps in pc 



17 
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0.4 
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Figure 7: The saturation curves, Dq, for various values of the gluon condensate. 
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Figure 8: The saturation curves, Dq, for various values of the vacuum saturation 
violation parameter, V^s- 
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range < c < 2. In Fig. ^ we see how deviations from vacuum saturation affect our 
results. We find reasonable agreement for < V^s < 2, suggesting that vacuum saturation 
may be violated by 100% with little effect. For the remaining OPE parameters, varying 
their values from to 200% was found to have no noticeable effect. However, the best 
results are obtained with the strong coupling taken to be close to the central value quoted 
in Table |l|. 



6.3 The Average Light Quark Mass 

Another way to measure whether the sum rules are saturated is by making a prediction 
for a physical quantity. As is an unphysical parameter, for 'sensible' values of sq this 
prediction should be effectively independent of over a wide range. From Eq. (0), we 
see that the theoretical side of our sum rules contains an overall factor of rn^il GeV^), 
we can thus make an estimate for the average light quark mass from 



m,(l GeV 



k\had 



(35) 



calc 



If the graph of this quantity against is found to be flat then we can say that the sum 




1.5 ,2 ^2.5 
m2 (GeV2) 

Figure 9: The quark mass curves for our two Coupling Schemes. The dashed line 
shows Coupling Scheme I and the solid line corresponds to Coupling Scheme II. 



19 



rules are saturated. Although the actual value of the quark mass extracted will depend 
on the constant /, its dependence will not, thus as a test of saturation Eq. (|35|) is 
independent of any uncertainties in the normalisation of the hadronic side. Once again 
we will use only the lowest sum rule, i.e. Tq, to estimate the quark mass. 



5.5 p 

5.4 - 

5.3 

f5.2 
S 

e"5.1 

5 



p. = (470 MeV)-l 
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(b) The Gluon Condensate 
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(c) Vacuum Saturation Violation Parameter 
Figure 10: The effect of various parameters on the quark-mass curves. 



In Fig. ^ we show the curve obtained from Eq. (|35D using the instanton parameters 
determined in Sect. for our two Coupling Schemes. We notice immediately that the 
curve is extremely flat, in particular for Coupling Scheme II. In both cases the variation 
is less than 1% in the range 0.6 GeV^ < < sq. This confirms that our sum rules are 
well satisfied even for quite low values of M^. The difference between the results of the 
two Coupling Schemes is ~ 1%. Thus we estimate mg(l GeV^) = 5.2 MeV. 
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Technique 


ruqil GeV^) (MeV) 


Chetyrkin et al. 
Prades 
Maltman & Kambor [52 
Dosch 53] 
APE 98 |54 
JLQCD 99 |55| 
QCDSF 99 

APE 99 |57 
CP-PACS 00 |58 
SESAM 98 |59| 

CP-PACS 00 ig 

QCDSF - UKQCD 01 || 


Pseudoscalar Sum Rule 
)) 

SR for (gg) + GMOR 
Quenched Lattice QCD 

» 

Unquenched Lattice QCD 

77 


5.7 ±1.2 
6.4 ±2.3 
5.6 ±0.8 
4.7-7.9 

6.8 ±0.6 
6.0 ±0.4 
6.3 ±0.3 

6.8 ±0.7 
6.2 ±0.2 

3.9 ±0.2 
4.9 ±0.3 
5.0 ±0.3 



Table 2: Recent determinations of mg(l GeV^). Note that Lattice QCD results 
are generally quoted as mg(4 GeV^) and have here been scaled up by a factor of 
1.42 as required by Eq. (|T^) with ^^(m^) = 0.334. 



We have performed a number of tests to determine how stable this estimate is to the 
variation of input parameters within the ranges found to give good saturation curves. The 
parameters pc, {agG"^) and V^s each yield an uncertainty in the quark mass of 0.1 MeV, 
see Fig. 0. Varying the remaining condensate values from zero to twice the value listed 
in Table [l| produces no noticeable effect on the quark mass obtained. 

The largest source of uncertainty in our calculation is the value of the phenomenological 
constant of proportionality, /, defined in Eqs. (27,28). This could only be determined to an 
accuracy of 20%, and so gives the dominant uncertainty in the final result of 0.5 MeV. Our 
final estimate is then mq{l GeV^) = 5.2 ± 0.6 MeV. In Table we list some other recent 
determinations of this quantity. We see that our estimate is consistent, though lower, 
than earlier sum rule determinations. Though not agreeing with the results of quenched 
lattice QCD, our value is reassuringly consistent with the early results of unquenched 
lattice calculations which are just beginning to appear in the literature. 

7 Discussion 

The application of QCD sum rules to the isoscalar-scalar current requires a realistic rep- 
resentation of the contributing resonances. None in this channel remotely approximates 
a 5— function. Fig. 2. Indeed, states like the a (or /o(400 — 1200)) are far too broad to 
be sensibly described by simple Breit-Wigner forms. Moreover they overlap with other 
resonances like the /o(980) and strongly coupled thresholds, like KK. It is quite natural 
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that the I = J = correlator is expressly related to the vrvr scattering amplitude and so 
the vacuum polarization can be described directly in terms of experimental information 
without the need for further (Breit-Wigner-like) approximations. We have shown that 
such Coupling Schemes, defined in Sect. are perfectly capable of saturating QCD fi- 
nite energy sum rules. In vrvr scattering the /o(980) appears as a sharp dip. Consequently, 
in our coupling schemes this state makes relatively little contribution to the correlator. 



in total contrast to the conclusions of Shakin and Wang ||6T[. That good saturation is 
possible strongly suggests then that the /o(980) does indeed make only a small contri- 
bution to the scalar-isoscalar sum rules, which in turn implies that it does not contain a 
large nn component in its wave function. This result is also at odds with the conclusion 



reached in |]13|, jl^. However, as explained above, when the same authors [|14| considered 
a more realistic approximation for the resonance shape, the parameters found were more 
consistent with the /o(400 — 1200) than the /o(980). From the present analysis, beyond 
saying that it is not strongly nn, we can make no comment on the structure of the /o(980). 

The saturation curves we obtain are remarkably flat and close to one from low values 
of M^, see Figs. 4-8. As the exponential in the weight function tends to suppress the 
integrand in the region s > M^, saturation at these values of means it is the lower 
energy portion of the experimental data that is most important in saturating the sum 
rules. Furthermore, the /o(1370) (with mass m and width F) will not contribute until 
> (m - F/2)2 ^ 1.4 GeV^. That the sum rules are fulfilled from 1 GeV^ implies 

it is the /o(400 — 1200) that plays the most important role in saturating the nn scalar- 
isoscalar sum rules. 

Our results again show that the Operator Product Expansion is not enough to describe 
the physics of the scalar mesons completely. We find saturation only when instantons 
are included. One potentially worrying point of this analysis is the need to tune the 
average instanton size away from its accepted value of pc ~ (600 MeV)"^ ~ 0.33 fm. to 
Pc ^ (470 - 490 MeV)-^ or 0.40 - 0.42 fm, which is about 25% larger. The validity of the 
Instanton Liquid Model relies on the mean separation of instantons being significantly 
larger than their average size. If this is not true then neighbouring instantons will overlap 
significantly and it becomes meaningless to talk about individual instantons. The assumed 
average instanton separation of 1 fm is only about 2.5 times the value of pc required by 
our sum rule, and so we are probably approaching this limit. Reassuringly though, lattice 
simulations of the topological features of the QCD vacuum seem to find average instanton 
sizes similar to the one we are forced to use (for a summary of lattice determinations of 
instanton parameters see p2|). 
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The above results follow from ratios of the finite energy sum rules defined in Eq. (33) 
and displayed in Figs. 4-8. Moreover, the absolute value of each is determined by the 



quark mass curves of Sect. |6.3| . These are also found to be fiat and this constancy sets in 
already below = 1 GeV^ Fig. 9. This is again evidence that our Coupling Schemes 
are capable of satisfying the sum rules and that the state most important for this is the 
/o(400 — 1200). The value we find for the average u and d quark mass is 

mg(l GeV^) = 5.2 ± 0.6 MeV . (36) 

The most conservative way to interpret this result is merely as a consistency check. A num- 
ber of assumptions are made in this analysis, e.g. the validity of the Coupling Schemes, 
the truncation of the OPE at dimension-6 and the perturbative contribution at 0{al), the 
form of the instanton contribution etc. That our final result is of the correct size reassures 
us that none of our assumptions are too far off the mark. Going further and taking it 
foregranted that the isoscalar-scalar correlator is indeed directly related to vrvr scattering, 
then the sum rules provide an independent determination of the average light quark mass, 
as given in Eq. (37). Our result is quite compatible with sum rule determinations in the 
pseudoscalar sector and with the early unquenched lattice estimates that are starting to 
appear in the literature. We can then be certain that the /o(400 — 1200) is the dominant 
contribution to the uu + dd current. 
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A Renormalisation Group Conventions 



The QCD (3 and 7-functions are defined by the the Renormahsation Group Equations 

i>0 



d/x2 



(37) 



and 



2 

37T = -mq-f{as) 



d/x2 



(3^ 



j>0 



where = a^/vr. 

Currently, the QCD /3 and 7-functions are known to four-loop order. For three active 



flavours, the coefficients of the (3 function within the MS scheme are 
9 



Pi 



P2 



3863 
^8? 



140 599 445^,^, 
P^ = ^7:777^ + — C(3) . 



4608 



32 



(39) 



The coefficients of the 7-function, also within the MS scheme, are ||64| , |65 

91 



7o = 1 , 7i 



24 



8885 5 



2977517 9295^,^, 135^,,, 125^,^, 
73 = -^T7^ - ^C(3) + — C(4) - — C(5) 



41472 



432 



16 



12 



(40) 
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